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Abstract 

We update our study of neutrino mass hierarchy determination using a high statistics reactor 
D e experiment in the light of the recent evidences of a relatively large non-zero value of #13 from 
the Daya Bay and RENO experiments. We find that there are noticeable modifications in the 
results, which allow a relaxation in the detector's characteristics, such as the energy resolution 



^vq I and exposure, required to obtain a significant sensitivity to, or to determine, the neutrino mass 

k£ hierarchy in such a reactor experiment. 

1 Introduction 



Determining the type of neutrino mass spectrum, which can be with normal or inverted ordering 
(NO or 10) or hierarchy (see, e.g., pQ), is one of the most pressing and challenging problems of 
future research in neutrino physics. The recently measured relatively large value of the angle #13 
of the Pontecorvo, Maki, Nakagawa, Sakata (PMNS) neutrino mixing matrix in the Daya Bay [2] 
and RENO [3] experiments opens up the possibility of the neutrino mass hierarchy determination 
in an experiment with reactor v e . This possibility was discussed first in [6] and later was further 
investigated in [91 [TO] (see also [11] ) . It is based on the observation that for cos 2#i2 7^ and 
sin #i3 7^ 0, 9\2 being the solar neutrino mixing angle (see, e.g., [1]), the probabilities of v e survival 
in the cases of NO (NH) and 10 (IH) spectra differ Q2]: P NH [p e -> v e ) / P m {v e -> v e ), and 
\P NH (p e — > i> e ) — P 1H (v e — > v e )\ oc sin 2 2^i3 cos 26*12. For sufficiently large | cos 2^12 1 and sin 2 #13 
and a baseline of several tens of kilometers, this difference in the v e oscillations leads, in principle, 
. to an observable difference in the deformations of the spectrum of e + [B] , produced in the inverse 
beta-decay reaction v e + p — > e + + n by which the reactor v e are detected. 

In the present Addendum we re-evaluate the potential of the reactor v e experiments for de- 
termination of the neutrino mass hierarchy using the Daya Bay and RENO data on #13. Such a 
re-evaluation is necessary since sin 2 #13 was measured with a relatively high precision in the Daya 
Bay and RENO experiments and found to have a relatively large value. We expect the latter to 
lead to less demanding, than previously estimated, characteristics of the v e detector, required for 
getting information about the type of the neutrino mass spectrum. 

We perform the analysis uing the methods descibed in detail in [TO]. We assume the experiment 
is performed with a KamLAND-like (see, e.g., [13J ) 10 kT detector (planned, e.g., within the project 
Hanohano [H]), located at L = 60 km from a reactor v e source, having a power of ~ 5 GW. As 
in [TO] (see also [7]), the threshold of the visible energy used is set to E vist h = 1.0 MeV. As is 
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well known, for the experimentally determined values of the solar and atmospheric neutrino mass 
squared differences, which we give below, the optimal baseline of the experiment of interest is 
approximately 60 km (see, e.g., [S1E]). We present results also for the shorter non-optimal baseline 
of L = 30 km. For the reactor angle #13, we use the results of the Daya Bay experiment |3j: 

sin 2 20i3 = 0.092 ±0.016 ±0.005, 0.04 < sin 2 26 13 < 0.14 , 3a . (1) 

In what concerns the other oscillation parameters which enter into the expressions for the reactor 
v e survival probabilities in the cases of NO and IO spectra, the solar and atmospheric neutrino mass 
squared differences, Am|, = Am^ and Am^ = Am^ = Am 2 2 , and the solar neutrino mixing angle, 
#12 , we use the values obtained in the global analysis of the neutrino oscillation data, including the 
data from the Daya Bay and RENO experiments, performed in [16J. It follows from the results 
obtained in [TO], in particular, that we have cos2#i2 > 0.28 at 3a . 

Since the sensitivity to the neutrino mass hierarchy of a reactor v e experiment depends critically 
on the value of the angle #13, we have redone our earlier analysis [10] taking into account the new 
data on sin #13, eq. ([I]), including the allowed 3a interval of values. In the following Section we 
present our updated analysis and results. 

2 Updated x 2 -Analysis of the sensitivity to the type of the neu- 
trino mass spectrum 

We perform a full x 2 analysis of the hierarchy sensitivity of a medium-baseline reactor u e 
experiment with a detector of the prototype of KamLAND, choosing the optimal baseline of 60 km 
unless otherwise stated. The hierarchy sensitivity is computed by simulating an "experimental" 
event spectrum for a fixed "true" hierarchy (we choose a normal hierarchy here, the difference being 
minimal if it is chosen to be the inverted one). A "theoretical" event spectrum is simulated with 
the other or "wrong" hierarchy. A standard Gaussian \ 2 is then obtained, which determines the 
confidence level at which the "wrong" hierarchy can be excluded. 

Our rigorous analysis involves optimizing the event binning to give the best sensitivity while 
being compatible with constraints of detector resolution, marginalizing over the neutrino parameters 
|Am 2 tm | and #13, and taking into account systematic and geo-neutrino uncertainties by the method 
of pulls (for further technical details of the analysis see [TO]). We have checked in [TO] that doing 
a marginalization over sin 2 #12 and Am^ over their present 3a ranges of sin 2 #12 = 0.26 — 0.36 
and Am 2 ! = 7.0 x 10 -5 eV 2 — 8.2 x 10~ 5 eV 2 [TO] does not significantly affect the results on 
hierarchy sensitivity, since they are relatively small variations. Hence, we have presented in [TO] 
the final results with the values of sin 2 #12 and Am^ fixed at their best-fits of sin 2 612 = 0.31 and 
Am^ = 7.6 x 10~ 5 eV 2 . We follow the same procedure here @. 

We present results for different values of #13, the detector exposure and the energy resolution. 
As was done in [TO] , a prior term is added to take into account information from other experiments 
on parameter uncertainties. We find that the uncertainties in the values of |Am 2 tm | and #13 play a 
crucial role in the sensitivity to the neutrino mass hierarchy, and hence the reduction in the allowed 
range of #13 as well as its increased value aid in hierarchy determination. We study the effect of the 
detector energy resolution, exposure, parameter marginalization and data binning using the new 
data on #13, eq. ([I]). 

4 We have made use also of the results found in [10] (see also [3 E]) that the inclusion of systematic and geo- 
neutrino uncertainties as well as of ~ 1% energy scale shrink/shift uncertainty (even if energy-dependent), has only 
a minimal effect on the neutrino mass hierarchy determination. 
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Figure 1: The hierarchy sensitivity (x 2 )stat as a function of the number of L/E bins, for fixed neutrino oscillation 
parameters, sin 2 2f?i3 =0.1 and detector's energy resolution of 3%, statistics of 200 kT GW yr, baseline of 60 km and 
different L/E binnings in the range L/E = 5 — 32 km/MeV. 



We consider the following error ranges for the two marginalized parameters: i) | Aro^ | is allowed 
to vary in the range 2.3 x 10" 3 - 2.6 x 10~ 3 eV 2 , and ii) sin 2#i3 is varied from 0.04 to 0.14, to be 
consistent with the 3a range found in the Daya Bay experiment. 

Figure Q] shows the behaviour of the x 2 sensitivity with an increase in the bin number for fixed 
neutrino parameters and an exposure of 200 kT GW yr, using sin 2 29\ 3 = 0.1, Am^JVff) = 
2.4 x 10" 3 eV 2 , A 31 (IH) = -A 31 (NH) + Am^ and a detector resolution of 3%, for different 
numbers of L/E bins in the range L/E = 5 — 32 km/MeV. The sensitivity is seen to improve 
dramatically with an improvement in the fineness of division, and the binning is optimized at 150 
L/E bins to derive the best possible sensitivity while being consistent with the detector resolution. 
For 150 (100) L/E bins, the bin width in energy in the case we are considering is about 68 (100) 
keV. 

Table 1 lists the values of the hierarchy sensitivity (x 2 )Sa™ f° r different values of #13 and the 
detector energy resolution, after a marginalization over the parameter ranges indicated above, for 
an exposure of 200 kT GW yr and a 150-bin analysis. The true #13 values are chosen within the 3a 
range allowed by the Daya Bay data. Prior experimental information regarding the other neutrino 
parameters is included in the analysis in the form of "priors", using the present la error ranges 
of the respective parameters: fj(|Am 2 tm |) = 5%x | Am 2 tm | truc and cr(sin 2 26> i3 ) = 0.02. Table 2 
gives the values of the hierarchy sensitivity [(x 2 ) 1 !^^^ for different values of 6*13 and the detector 
energy resolution with a parameter marginalization including priors, for the same values of detector 
exposure and event binning. The slight improvement in the results with the inclusion of priors is 
enhanced if a lower prospective la error of <r(sin 2 28\ 3 ) = 0.01 is considered. As recent reports 
from Daya Bay and RENO have shown, such an improvement in the precision of #13 is not far out 
of reach of present experiments. Moreover, a combined analysis of the global data on the angle #13 
performed in [15] already yields <7(sin 2 26> 13 ) = 0.013. 

In Table 3, we list the values of the hierarchy sensitivity [(x 2 )X?] f° r sin2 2 ^i3 Ue = an d 0.1, 
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Table 1: Values of (x 2 )S"t marginalized over the parameters 6 13 and |Am§i|, for lAm 2 !!' 1110 = 2.4 x 10~ 3 eV 2 , 
crfJAmli]) = 5%x|Am3 1 | truc , cr(sin 2 29 13 ) = 0.02, three values of sin 2 26\ T 3 ue and three values of the detector energy 
resolution. The detector exposure used is 200 kT GW yr. The baseline is set to 60 km. The values of (x 2 )^!™ are 
obtained in an analysis using 150 L/E bins in the range 5-32 km/MeV. 



[(X )stat]prior 


Energy resolution 






sin 2 26\f c 


2% 


3% 


4% 


0.07 


6.37 


5.15 


3.90 


0.1 


13.17 


10.58 


8.05 


0.12 


19.10 


15.26 


11.60 



Table 2: The same as in Tabled] but for cr(sin 2 20i 3 ) = 0.01. 



for 3 different values of the detector resolution and a scaling in the detector exposure. These results 
show the strong dependence of the sensitivity on the detector exposure. For example, a hierarchy 
sensitivity of nearly 3a may be possible even for sin 2 26^ lc = 0.07 and an energy resolution of 
4%, with an exposure of 400 kT GW yr, and this would improve further with a higher detector 
mass/power. 

To highlight the improved sensitivities possible even for smaller detector exposures when #13 
is close to the present best-fit value, we present in Table 4 the hierarchy sensitivity [(x 2 )™at] f° r 
sin 2 26^ r 3 ue = 0.1 and 0.12 for lower detector exposures 100 kT GW yr and 150 kT GW yr with 3 
different values of the detector's energy resolution. We note that even with an energy resolution of 
4%, a 2(j sensitivity is achievable with a relatively low exposure of 100 kT GW yr for the indicated 
values of 6*13. With a better energy resolution of 2%, the sensitivity can go up to 3a or even to a 
higher value. 

In Table 5 we list the values of hierarchy sensitivity obtained for two detector exposures and 
three values of detector resolution when the baseline is chosen to be 30 km instead of 60 km. This 
table shows that the sensitivities decrease for the indicated shorter baseline, i.e., when the baseline 
deviates significantly from the optimal one of 50 - 60 km. For example, with a baseline of 30 km, 
a resolution of 2% and an exposure of 200 kT GW yr would be required for a hierarchy sensitivity 
of 3a if sin 2 26^ e = 0.1, while with a baseline of 60 km similar sensitivity is achievable with an 
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Table 3: The same as in Table[T] but for three values of the detector exposure and sin 2 2# 1 3 UC = 0.07; 0.1. 
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Table 4: Values of (x 2 )^ta" marginalized over the parameters Q13 and |Am§i| for lower detector exposures (in kT 
GW yr), sin 2 26 , 13 uc = 0.1 and 0.12, for three values of the detector's energy resolution and a baseline of 60 km. The 
results are obtained in an analysis using 150 L/E bins in the range 5-32 km/MeV. 



exposure of 150 kT GW yr with the same detector resolution [3 



3 Conclusions 

We find that the data on the parameter #13 from Daya Bay experiment allow us to get information 
or determine the neutrino mass hierarchy with a greater efficiency, than was previously estimated, 
using a reactor u e experiment: the stringent requirements of the detector's energy resolution and 
exposure obtained in the previous studies can be relaxed significantly. Since hierarchy sensitivity 
depends strongly on the the true value of the energy resolution and the exposure, a relatively 
large value of sin 2 2#^ 3 ue close to the Daya Bay best fit of 0.092 makes it easier to achieve hierarchy 
determination using lower detector exposures and less demanding energy resolution. 

For example, (x 2 )Sa? f° r the "wrong" hierarchy improves from 3.5 (1.8<t sensitivity) for sin 2 2#i3 UC 
0.05 (close to the Daya Bay 3a lower limit), an energy resolution of 2% and a detector exposure of 
200 kT GW yr, to 12.9 (a 3.6cr determination) for sin 2 26>f 3 ue = 0.10 (close to the Daya Bay best fit) 
for the same values of the resolution and exposure. With this value of sin 2 2#^ 3 uc , even an energy 
resolution of 4% can give a sensitivity of nearly 3a. 

To summarise, for the values of #13 from the interval allowed at 3<r by the Daya Bay data, a 
significant hierarchy sensitivity is possible even with a detector energy resolution of a ~ 4% and 

5 The optimal baseline for hierarchy sensitivity lies in the region of maximization of the effect of 
the phase Am| 1 L/2i5 in the expression for the v e survival probability. With the present error range of 
Am 2 !, and the peak of the reactor v e event rate spectrum at 3.6 MeV, this gives an optimal baseline 
range of 55 to 64 km. Hence, the hierarchy sensitivity becomes worse for baselines significantly shorter 
than the indicated range. 
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4.90 


3.80 
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8.79 


6.50 


5.05 


12.81 
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7.35 



Table 5: Values of (x 2 )™at marginalized over the parameters #13 and lAm^l for two values of detector exposures 
(in kT GW yr), for sin 2 2^i3 UC = 0.1 and 0.12, three values of the detector's energy resolution and a baseline of 30 
km. The results are obtained in an analysis using 150 L/E bins in the range 5-32 km/MeV. 



an exposure of 200 kT GW yr. For sin 2 26^3 UC = 0.10 (0.12) and an energy resolution of a ~ 2%, 
a 3a sensitivity to the neutrino mass hierarchy can be achieved with an exposure of 150 (100) kT 
GW yr. The indicated requirements on the detector specifications make the discussed reactor v e 
experiment more feasible than the previous analyses have suggested. 
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